Finite temperature damping of collective modes of a BCS-BEC crossover superfluid 
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A new mechanism is proposed to explain the puzzling damping of collective excitations, which was 
recently observed in the experiments of strongly interacting Fermi gases below the superfluid critical 
temperature on the fermionic (BCS) side of Feshbach resonance. Sound velocity, superfluid density 
and damping rate are calculated with effective field theory. We find that a dominant damping 
process is due to the interaction between superfluid phonons and thermally excited fermionic quasi- 
particles, in contrast to the previously proposed pair-breaking mechanism. Results from our effective 
model are compared quantitatively with recent experimental findings, showing a good agreement. 
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I. INTRODUCTION 

In recent years, experimental studies on ultracold 
Fermi gases in Bose-Einstein condensation (BEC) to 
Bardcen-Cooper-SchriefFer (BCS) crossover [l|-yl have 
shown many interesting results [J-|7|. By varying exter- 
nal magnetic field to effectively tune the scattering length 
through Feshbach Resonance, or changing the tempera- 
ture, collective behaviors, such as oscillation frequency, 
sound velocity and damping rate, have been extensively 
investigated |8l-[l3|. Recently, several experiments show 
a puzzling damping of collective excitations occurs in su- 
perfluid regime where the system is at a finite temper- 
ature and on the (fermion) BCS side of the crossover 

In the theoretical front, the properties of collective ex- 
citations in the crossover have been investigated by a 
variety of established methods 3, in cluding work that 
used the hydrodynamic theory |13l - [l5| , the Gaussian fiuc- 
tuation approximation from a microscopic Lagrangian 
p^| - [l9j . and fermion-boson coupled model |20l - [2^ . The 
Boltzmann equation |23l - l26l | and semiclassical dynam- 
ics [23, [23I were also used to explain the frequency and 
damping rate of oscillation modes. 

Effective field theory has also been used in some pre- 
vious work to study ultracold systems at zero tempera- 
ture, including expanding the goldstone fields to higher 
orders, or low energy expansion J29l - l3l| . In this pa- 
per, we use effective field theory to study the collective 
modes of the superfluid state of an unpolarized strongly 
interacting cold Fermi gas at flnite temperatures near 
the Feshbach resonance regime. Our effective theory is 
constructed based on a single-channel resonance model, 
which is known to be adequate for a broad Feshbach res- 
onance. Furthermore, we focus on the BCS side of the 
crossover, where puzzling damping phenomena have been 
observed and received much attention [^ [III, [l2| • In this 
regime, the BEC molecular effect is negligible, and hence 
the single-channel calculation which we adopt is simple 
and valid. Our main results are as follows. First, our the- 
ory shows that a dominant damping process is due to in- 



teraction between superfluid phonons and thermally ex- 
cited fermionic quasi-particles at finite temperatures, in 
contrast to the previously proposed pair-breaking mecha- 
nism. Such a process is related to the well-known Landau 
damping, which was previously discussed in a semiclas- 
sical approach by taking the S — > limit of the time- 
dependent Hartree-Fock-Bogoliubov equations [23, [28| . 
Second, we calculate the damping rate of collective exci- 
tations for the physical systems and find a good compar- 
ison with the recent experimental findings. Throughout 
the paper, the Planck constant h and Boltzmann con- 
stant ks are set unity (ft = fc^ = 1 in units). 



II. MEAN FIELD THEORY 

We consider an unpolarized ultracold fermionic system 
with two spin species (spin f and ],) of equal mass. The 
partition function of the system in the grand canonical 
ensemble can be described in the path integral formalism 



Z = I D{r,,^,)eM-S^) 



(1) 



where S*^ = J dxCf{x) is the action, ip^r is the fcrnfionic 
field for spin a =t,i- a; is a four vector x — (x, r), 
in which x is the spatial coordinate and r is imaginary 
time in the range < r < /3, /3 = 1/T and T is the 
temperature of the system. The following Lagrangian is 
used to describe the system: 



^4> = V'ct 5r 



2ma 



- n]i'a+ g'ip^'^Pl^Pi1p^, (2) 



where rUa is the mass of the fermion, g is the interac- 
tion (negative on the BCS side), and fi is the chemical 
potential. The summation over a is implicit. 

As usual, we use a Hubbard- Stratonovich transforma- 
tion by introducing an auxiliary complex bosonic field 
A(a;) to eliminate the quartic term in Eq. ([2]) and get 



Z ^ 



y Z?(C, ^^)^(A*, A) exp(-5^,A 



(3) 



where the lagrangian is now 



I^^A = V'ff c^T 



/iU-^ + C^IV-IA + cc)- 
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(4) 

In this transformation, the (auxihary) field A(x) acquires 
exactly the same expectation value as the pair field at the 
saddle point: 

(A) = (V-i^t) ■ 



By self consistently solving Eq. (fTTj) . we can get Aq 
and jjL as function of 1/kpa and T/Tp in mean field level, 
where kp and Tp arc the Fermi momentum and temper- 
ature of the non-interacting Fermi gas in free space with 
the same number density n. In the following sections, we 
will assume the state of the system is already superfluid 
and will use the mean field results, such as the gap pa- 
rameter, as input to calculate the physical quantities of 
excitations. 



The fermionie field is now in quadratic form and can be 
integrated out to get an effective action for the field A 



Z = J D{A*,A)cxp{~Sa), 



where 



and 



Sa 



d4^f ^+trlnG-i 



G- 




(5) 



(6) 



(7) 



In the mean field approximation, we seek for a spatially 
uniform saddle point solution Aq of 5a, together with 
the requirement of fixing the number density n. These 
two conditions determine the mean field equations in the 
crossover 



5Sa 

SAo 

dlnZ 



0, 



(8) 



In momentum space we get 

— h >^ tanh 
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= 0, 



1-iitanhf^ 
Ek V 2 



(9) 



where E^ = VlT+A^ and S,k = ^ - IJ.- 



We use the usual regularization procedure, 

1 

2^ 



rUa 1 Y^ 1 

4^ ~ g(A) "'' ^ '^ 

^^ ' k <A 



(10) 



where Ck = |kp/2ma and a is the effective scattering 
length. Taking A — > oo, Eq. (|9]) can be expressed in 
terms of a as 



47ra ^ 



E 



^-^tanhf^ 
2efc 2Ek \ 2 



1-i^tanhf^ 
Ek V 2 



(11) 



III. QUANTUM FLUCTUATIONS 

To consider the quantum fluctuations, we no longer 
treat the order parameter A as a homogeneous constant 
in Eq. Q, but write it as 



A(a;)== Ao(l + A(a;))e'2v^(^). 



(12) 



A is the amplitude fluctuation and (/3 is the phase fluctu- 
ation around the saddle point solution Aq, and both of 
them are real. Then one can apply a local U(l) gauge 
transformation to a new gauge where the order parame- 
ter is real everywhere in space 



A = Ae*2v(^) , Va (x) = i^a (a;)e'^(^) , 



(13) 



and A = Ao(l -I- A) is real now [S^l- In the new gauge, 
the Lagrangrian density Eq. (|4]) becomes 



C 



■Tp.'^.X 



rAdr 



2?Tin 



fl)'lpa + (AoV'tV'l + C.C.) 



2ma 

+{A(i\^i'l + C.C.) 



Al\^ 



(14) 



where J, = J^x) - -^^i^A^^^) - (VC)V'^] is the 
fermion current field. 

We now integrate out the fermionie field which is 
gapped. In momentum space, for the fermionie field 
ijj, we use k — (k, iojm), where aj,„ = {2m + l)7r//3 is 
the fermionie Matubara frequecy and m is integer. For 
the bosonie field A and ip, we use q — (q, icj„), and 
Ldn = 2n7r//3, where n is integer. 

Introducing two component spinor in momentum space 
^k ~ i''Pk■\■^''P-kl)^ '^6 can rewrite the action in mo- 
mentum space as 5" = J2k k ^l G^^{ki,k2)'ik2J where 
G~^ = Gq^ + xi + X2 is 2 X 2 matrix. Go, Xi, and X2 
are given by 

ri n 1 \ ^kiM ( i^mi+iki Aq 

Go(fcl,fc2)=^,2 ^P2 I A„ 



^mi + K 



i^mi - Cfci 



!i I, \ >r- q- (ki -ka -q) / 1 

X2(kl,k2) = - 2_^ ^_ fl-,^ fqfki-k2-q 



2mapV 



Xl(fcl,fc2)-(^^21(fc^^fc^) Xf(fcl,fc2) 



-1 
(15) 



with 

X"(fcl,fc2) 



2 _1,2M 



zikf - ki 






i^rni ^-^7712 



V^ 2mV^ 






xlHkuh) = xf{h,k2) 
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/^ 



Afei-fe2- 



Keeping the A and (/? field to quadratic order [33 



Z?(^* vl/)e~'^*'*''^ sa e'''[*^''-^^~5'^"-^i'~^''-^^i = e 



we get the effective action 

i^ ^ Otp -\- ox + t^ipjA? 

where 

'S'tp = ^ '/'9'/'-g[/i(9)w^ + lUiq) 



2m„ 



q 
S^.x = Y][/2(g)A5(^_5(-w„ q • k) 



-f2{q)X-qfq{~^n H Q ' k)] , 



and m 



1 ■^-^ WniCLim-i-n — £,k^k+q + Ag 






g /3F ^ (a.?„ + El ) (..2^+„ + El^^ ) ' 

J_^ fc^ COS^ 9 ~LL>m(^rn+n + CfcCfc+g + '^0 



(20) 



The next step is to integrate out the amphtudc fluctu- 
ation field A and keep only the phase field. Define 



Cl(g)-/l + /|//3, 
C2{q) = .U + 



2ma' 



(21) 



and the effective action of phase field ip is obtained as 

5^=^^,^_,[Ci(g)c^2^C2(<z)q2]. (22) 



IV. SUPERFLUID DENSITY AND SOUND 
VELOCITY 



Some physical quantities such as superfluid density and 
sound velocity can be evaluated by our theory. To get a 
full expression for the superfluid density, one can apply 
a Galilean boost on the order parameter in Eq. p^ . and 

(16) get the superfluid density by the shift of free energy in 
the superfluid hydrodynamic model with normal and su- 
perfluid components [13, [la Sj • The total normal den- 
sity includes bosonic fluctuation part and non-condensed 
fermionic Bogoliubovquasi- particle part. From the cal- 

-s' culations of Ref. [IsllajlSalj the main contribution to the 

' normal density on the BEG side is the bosonic fluctua- 

(17) tion part, while on the BGS side, the main contribution 
is from the non-condensed fermionic quasi-particle part. 
Since our main focus is the damping on the BGS side, we 

(18) will neglect the part of bosonic fluctuations in the normal 
density in the following calculation. 

To the zcroth order of g, or in the long wave length 
limit g — > 0, all coefficients /i,/2,/3, and f^ are real, 
in terms of which superfiuid density and sound velocity 
can be expressed. Since the imaginary part is zero, there 
is no damping for the collective modes, to the lowest 
order. C[ = limq_i.o Ci is related to the density of states, 
and C2 = limQ^o C2 — ns/2ma is related to superfluid 
density Us [36[ J37| . 

Take the limit of g -^ in the coefficients Eq. (|20p . 
and define /( = limg_>g /i and so on. After carrying out 

(19) the Matsubara summation for iujm, we get [38| 
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/2 = # Et/^^'^-K - 1)^ + (1 - 2n.) J3] ' 



V 



f'^-^ll'^PMnk-l)^ 
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2El 
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V mi 



(23) 



where 



nk 



o0Ek 



1 



is the quasi-particle Fcrmi-Dirac distribution. Eq. (|22 
becomes 



with 



5 = E^«^-?[^i^" + ^2q' 



2ma 



(24) 



(25) 




1/kpa 



FIG. 1. (Color online) (a) Superfluid sound velocity compared 
with Fermi velocity Vs/vp in superfluid regime, (b) Superfluid 
density compared with the total density Us/n. 



LUq + iO"*", the zero 



After applying a Wick rotation iw„ 

of the action, Eq. ([24| . gives the spectrum ujq = Ws|q|, 
where Vs is the sound velocity. Thus, we get the following 
physical quantities. 



n 



2ma 



v' = 



2ma 
/4 + n/2ma 



At zero temperature, Eq. ((23)) reduces to 



(26) 



n 



T=0 



J2.T=0 ■" 



V ^2i?f 

V ^ 2EV 



. _A4 1 



V 
0. 



(27) 



Perform the summation for /"s over k as integration in 
the thermodynamic limit, and define the following inte- 
gral quantities 



J2 = 



J4 = 



h 



dk— 

El' 



dk 



After some calculations, the sound velocity at T 

1 J2J4A2 



J)2 



3m2j|A2 + J|' 



(28) 
Ois 

(29) 



which reproduces the result of Ref. [39[ , albeit that was 
derived in a different method. 

Fig. [T] shows the numerical plot of superfluid density 
and superfluid sound velocity with varying temperature 
and scattering length in the BCS regime. 
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FIG. 2. (Color online) Two separate damping channels. The 
blue curves are the fermionic quasi-particle spectrum. At 
T = 0, the lower branch is fully occupied and the higher 
branch is empty. The green arrow is the incident phonon. 
The black arrow indicates the excitation of a fermionic quasi- 
particle. The red straight line is the linear spectrum of the 
phonon, where the slope is the sound velocity. The two black 
dashed lines are the changes of momentum and energy of 
fermionic quasi-particle, as indicated, (a) A phonon scatters 
an existing (by finite temperature) Fermi quasi-particle to a 
different state in the same energy branch. Such a process is 
known as Landau damping. This process can also happen in 
the lower branch at finite temperature, which is not shown 
here, (b) A phonon creates a quasi-particle and -hole pair 
across the lower and upper bands, equivalent to Cooper pair 
breaking in the representation of original particles. 



DAMPING 



A. Approximate formulas of the damping rate 



We start from the effective action, Eq. (|22l) . for the 
superfluid phase field. Keeping q small but finite when 
evaluating the coeffiencets Ci and C2 in Eq. (PT|) , imagi- 
nary terms appear in Eq. (j20p , which corresponds to the 
damping of the collective modes. To get exact dispersion 
relation including damping for the superfluid phonons, 
one needs to sclf-consistently solve for poles of the effec- 
tive action Eq. ^F^ of (/?. However, in the regime where 
the damping is small, we can simplify the calculation. 

The dispersion relationship is determined by setting 
the quadratic field term Eq. (P^ to be zero. To evaluate 
the damping rate we apply a Wick rotation iw„ = loq — i^ 
as follows, 



-Ci(9)(c^o-*7)' + C2(gV 
q= (q, iun -^ wo + iO+), 



0, 



(30) 



where wg is the oscillation frequency, or the real part 
of the phonon mode, and 7 is the damping rate. Al- 
though in principle, we should also keep 7 when apply- 
ing the Wick rotation to the w„'s in coefficients /i, /2, /a, 
and /4, for small damping it is sufficient to just apply 
iujn = wq + *0"'" for the coefficients to get the lowest order 
results. This approximation is similar to that adopted in 
the calculation for self-energy, for example, Ref. [4^ . 

For the numerators in Eq. ([20]) . we can still set g = 0. 
The reason to do this is that the numerators arc all real, 
and hence the small q expansion gives only high order 



corrections to the zeroth order contribution. If we want 
only the leading order in the long wavelength expansion 
of the damping rate, the high order corrections in the 
numerators can be ignored, as we have checked. 

Also, we can still use the results Eq. p3)) obtained 
from last section for the real parts of Ci, C2, /i, /2, /s, 
and /4, because wc arc not interested in the high order 
corrections to the real parts of coefficients, which just 
give higher order correction to the damping rate. 

Adopting the above approximation, we can write Ci = 
C[ + iC[', C2=C^ + iC'^, A = /{ + if I and so on, to find 
the poles of the effective action 

{C[ + iC'Oiuo - iif - [C^ + iC'^W = • (31) 

Solving this equation gives the damping rate 

7 _ 1 (C'{ 



Wo 



C'^ 



2\C[ C'^ 



(32) 



From Eq. (j2ip . assuming that all the imaginary parts are 
small, we have 



^1 — Jl + ^77 •'2 ^77 / 3 ' 



^2 — Ji ■ 



(33) 



After some calculation, we find that the imaginary 
parts of the coefficients /"'s can be separated into two 
channels a and fe, where in channel a 



vY 



Ei 



f'la = -^ E l^-^^t^" - (^'^■+« - ^^)K"fc - ^k+'i) ' 



fll 



fi.a 



V 

ttA, 



4 I 

- E T^^l-^o - {Ek+q - Ek)\{nk - rik+g) , 



^ V^^ 



V 2 '^['^0 ^ (^fc+9 ^ Ek)]{nk - Uk+q) , 

(34) 



and in channel b 



/l,6 



V ^^ 2E 



^6{ljo - Ek - Ek+q){l -rik- nu+q) 



l2,b 



k *; 

A2^J A2 



V ^ 2E\ 



-— 2(5(a;o - Ek- Ek+q){l - rik - rik 



+q) ; 



fib ^ E^O ^ 7^)'^('^0 - Ek- Ek+q) , 



V ^^"2 2Ei 

X(l-nfc -Uk+q), 

(35) 

For the expressions above, we would like to remind the 
reader of the definitions of 

1 



Q + A^ and Uk 



The physical meaning of channel a and 6, as defined 
in the above expressions, is illustrated in Fig. [51 Chan- 
nel a (Fig. [2]^a)) is the Landau damping. A (fermionic) 
quasi-particle absorbs or emits a superfluid phonon and 
becomes another quasi-particle state at a different mo- 
mentum within the same energy band. In channel b 
(Fig. [2][b)), a superfluid phonon excites a quasi-particle 
from the lower band to the upper band, creating a 
particle-hole pair excitation in the quasi-particle eigen- 
vector basis, which corresponds to a Cooper pair break- 
ing process in terms of original fermions. Both channels 
cause the superfluid phonons (collective exicitations of 
the superfluid state) to decay. 

Before applying our theory to real experimental sys- 
tems, let us discuss the properties of the damping in the 
wq — > limit (with wq/IqI = Vs flxcd) at low tempera- 
ture, where channel b vanishes. To satisfy the 6 function 
in Eq. (p4l). we need to have 



I cos 6*1 



Ekmvs 



ikk 



<1, 



(36) 



which is difficult to solve analytically. 

We can still get some analytical properties if we com- 
bine the numerical and analytical analysis. At l/kpci = 
—0.45 and T/Tp < 0.1, numerical analysis shows that 
k < kfj, does not contribute to the integral. For k > k^, 



we can substitute fc = /c^ in Eq. ((36)) . except for ^k, to 
get the range of ^k ■ We then apply a one-step iteration 
to get a more precise range for S^k- To satisfy Eq. ((36|) . 
we need to have 



a> 



VefTA^rm^ 



io = 



The next step is to write 



(37) 



{Ek - Ek + q) 



Uk - Tlk+q 

duk 

— • ^£.fc - r.k + q) 

l3uJonk{l - Uk) 



(38) 



The second line comes from ujq <^ T and the fourth line 
comes from T ^ A. Also, numerical analysis shows that 
the contribution from C" is much smaller than that from 
C2 in the above regime. In addition, at low tempera- 
ture, C2 ~ n/2ma since the normal density is negligible. 
Taking all the above into account, we get 



1_ 

Wo 



3ttv'! 1 



4 ^'z r Je,>^/^Ta^ 



dEk$e-^^'' 



(39) 



al3Ek 



The integral in Eq. ([5U)) needs to be evaluated numeri- 
cally in general cases. In the extremely low temperature. 



we can approximately treat E^/S,^. in Eq. p9p as a con- 
stant and substitute in S^j, = ^i , which gives 



3^«3(e2 + A2)2 ^A!T^ 



Wo 



"/ 



€t 



(40) 



An interesting feature of damping is that in Eq. (p9)) 
and (|40p. 7/1^0 is independent of wq when wg <C T. 
Therefore, we should observe the damping even in the 
Wo — > limit in finite temperature. A physical under- 
standing of this feature is that the damping from channel 
a is the coupling between superfluid sound (phonon) and 
thermally excited Fermi quasi-particles, so that 



7 ex dkf{k,A,fi){nk~nk+q) 



dK, 



ex Wo / dkf{k,A,iJ,) 



(41) 



where f{k,A,fi) is an arbitrary function independent of 
Wo. This provides an experimental method to verify our 
theory, which will be discussed later. 

In the following section, we will apply our theory to 
real experimental systems. We will focus on the regime 
of 7/wo ^ 1 on the BCS side of the Feshbach resonance, 
where our approximations are quantitatively controlled. 



B. Comparison with experiments 

The above theory can be applied to real physical sys- 
tems [41|. We apply our method to the experimental 
configuration in Ref. |ll| and calculate the damping rate 
as function of both temperature and scattering length. 
The number of particles is iV = 4 x 10^. The trapping 
frequencies are uj^ = 2tt x 830 Hz, Wj, = 27r x 415 Hz and 
w^ = 27r X 22 Hz. Thus^^w = (wj^Wj^w^)!/^ = 27r x 196 Hz. 
Reading from Ref. [il|, the oscillation (phonon) fre- 
quency is Wo ss 27r X 940 Hz sa 4.8 w. According to 
io/Tp = 1/(3A^)^/'^, we get the ratio between the phonon 
energy and Fermi temperature 



^=0.045. 
Tf 



(42) 



From Wo = w.s|q| and Ep = ^vpkp, we get 



|q| _ 1 Up wo_ 

kp 2 Vs Ep 

Since in most regime vp/{2vs) ~ 0(1) as shown in Fig.[Tl 
using Eq. (|^^ we conclude that |q| is also much smaller 
than the Fermi momentum kp. Thus, our model, which 
requires q to be small, can be applied to this physical 
system. 

Fig.inishows the numerical results of damping rate from 
channel a. It shows that the damping rate increases for 
higher temperature and smaller scattering amplitude \a\. 
In the superfluid regime, when T is close to T^, the damp- 
ing rate 7/wo becomes big. When solving the equations 




-0.12 



FIG. 3. (Color online) Damping rate of collective excitations 
in superfluid regime showing dependence on T/Tp and l/kpa. 
When T is close to Tc, the results from our approximation 
are no longer reliable since the damping rate 7 is already 
very large compared with cjq. This regime is indicated by 
the plateau. Also, when the system is no longer superfluid 
(T > Tc), our effective field theory does not apply, and the 
damping rate is then not plotted there. 



for the damping rate in the preceding section, we as- 
sumed a small damping 7 compared with the oscillation 
frequency wo, i.e., kept only the leading order in the per- 
turbative expansion of 7/wo. In the regime of big 7/wo 
where a significant correction is expected, our damping 
formula is no longer reliable. Thus, we use a plateau to 
indicate the regime of 7/wo > 0.5. Also, in the normal 
Fermi Liquid regime T > Tc where our theory no longer 
applies, we keep the damping plot open with no data 
points shown, and will discuss this regime later. 

We did not include channel b when plotting Fig. [31 The 
reason is that we found that under the experimental con- 
dition jll| , the contribution to damping due to channel b 
is much smaller than due to channel a. To investigate the 
features of channel a and channel b in more details, let us 
start from the relationship Wg — Ws|q|- We focus in the 
regime of long wavelength and low temperature such that 
kp > |q|, and T is not close to Tc. Thus A(r) is of the 
same order as Ep. In channel a, the first requirement is 
finite temperature, so that the upper (lower) band is pop- 
ulated by quasi-particles (quasi-holes) due to thermaliza- 
tion. This fact is enforced by the factor {nk — rik+q) in 
Eq. (|34[) . The second requirement is energy conservation. 
The energy change of the fcrmionic quasi-particle after 
scattering with a phonon is as follows, 

SE^E,,,-E,^^^<^\^M, (43) 

Ek nia Ek nia 

LOq « SE is needed to satisfy the energy conservation. 
Therefore, we have 



^'.|q| < 



6 |k||q| 



(44) 



Ek ma 

If |k| is too far away from kp, both the finite- 
temperature occupation number and the density of states 



arc greatly suppressed. Equivalently, the most effective 
scattering of phonons is from quasi-particle states around 
the upper (lower) band's minimum (maximum). In addi- 
tion, we have learned in Fig.[l]that Vs does not exceed the 
order of vp- As long as |k| is of the order kp or smaller, 
in principle, phonons of both large and small Vg may ex- 
cite fermionic quasi-particles in channel a. However, the 
condition Eq. (j44p suggests that small Vs is much more 
favored in channel a. 

In channel b, the pair breaking process is allowed by the 
condition of occupation number at T = (rik = nu+q = 
in Eq. ([35]) ). However, in this case. 



5E = Ek+g + Ek>2A, 



(45) 



where the minimal value of 6E takes place for k around 
kp. Again, by energy conservation, ujq « 5E is required. 
Subsequently, we find the minimal condition required of 
the sound velocity Vg, 



^ A kp 

Ek |q| 



(46) 



As long as A is of the same order of Ep , and given that 
fci?/|q| is quite large, the condition Eq. (l46l) in turns re- 
quires that the sound velocity Vs be much larger than 
the Fermi velocity vp. In the superfluid phonon case, 
the condition kp 3> |q| is satisfied, and when the tem- 
perature is still far away from Tc, channel b is prohibited, 
since Vg < vp as shown in Fig. [1] Another way to un- 
derstand that the damping channel b is suppressed is to 
directly use energy conservation. Since LOq <C 2A, chan- 
nel b can not happen. 

At finite temperature, channel a becomes possible since 
some fermionic quasi-particles and holes are thermally 
created in the upper and lower band, respectively. They 
can scatter with superfluid phonons to cause decay. How- 
ever, channel b is still greatly suppressed as long as 
Wq < 2A. 

In the experiment. Since tuo/Ep « 0.04, channel b can 
happen only when when A{T)/Ep « 0.02. From mean 
field analysis, near 1/kpa w —0.5, to satisfy A{T)/ Ep k, 
0.02, we need to have T /Tc ~ 99.7%, which is very near 
the phase transition. However, according to Ref. |ll| . 
the damping happens at T sa 0.6Tc near 1/kpa ~ —0.5, 
which is too low to let channel b happen from our calcula- 
tion. Even if one uses some theories including quantum 
fiuctuation [11], at T sa 0.6Tc, A{T)/Ep is still much 
larger than 0.02. Thus the damping peak should not cor- 
respond to the pair breaking channel b. As T getting 
closer to Tc, channel a will be more and more enhanced 
because of more and more thermally excited fermionic 
quasi-particles. When A{T)/Ep sa 0.02, channel b also 
happens, while the damping from channel a is already 
very large. It is not clear which channel dominates be- 
cause T near Tc is outside the valid regime of our low en- 
ergy effective field theory. It remains to be a challenge to 
formulate a quantum theory beyond the classical Boltz- 
mann equation. Thus, our calculation just considered 



the contribution from channel a. We also double checked 
channel b by numerical method and confirmed that the 
contribution from channel b is zero for most regime. 

There is another experimental evidence showing that 
why channel b does not dominate. In the experiments 
varying the magnetic field d, [l2|, if the pair breaking 
mechanism had dominated, one should also have ob- 
served very sharp peaks in these experiments. However, 
the damping rate changes relatively smoothly d Il3 |. 
which means channel a should be the reason for the 
smoothly increasing damping. Therefore, we expect the 
damping observed by changing the temperature |ll| to 
be due to channel a too, since increasing the magnetic 
field at a finite temperature has the same effect as in- 
creasing the temperature at fixing magnetic field, both 
just reducing the gap. Nevertheless, the pair breaking 
channel b is also possible to happen near the phase tran- 
sition, but it is not necessarily dominating in contrary to 
what has been suggested. 

The mechanism discussed here is different from the 
conventional acoustic attenuation process in solid-state 
superconductors, where channel b overwhelms channel a 



and 



2A. For the conventional case, Vs is big, so 



that a relatively large phonon energy ujq corresponds to 
a very small momentum q. Channel a is suppressed as 
Vs is too big to satisfy the energy conservation condi- 
tion (jm. At the same time, a channel b pair break- 
ing process (Fig. [DJb)), associated with small momen- 
tum but large energy transfer from acoustic phonons to 
Fermi quasi-particles, may happen. The original fermion 
pair breaking process is the creation of a pair of parti- 
cle and hole in the upper and lower branches of quasi- 
particle energy spectrum, respectively. A small |q| en- 
sures that fermionic quasi-particles are created at the 
band cxtrema, which is known to result in a peak of 
damping rate, 7 oc \/\/uj — 2A, due to the singularity 
of density of states. Therefore, the damping for large 
Vg in traditional case is due to the pair breaking process 
(channel 6). 

We did not consider the fact that in the experiment, 
the fermionic gases are trapped and inhomogeneous in 
space. The above calculation is effective for gas at the 
center of the trap. On the edge of the trapped gas, 
by local density approximation, the effective density and 
Fermi energy is smaller than that of center, which means 
ujo/Ep is larger and channel b may happen in lower tem- 
perature (but still not very low yet). However, channel 
a also happens on the edge since the above analysis still 
works for lower gas density on the edge, and the effect of 
channel b is just to increase the damping rate, not giving 
a sudden peak. 

In the regime T > T^ where the system is no longer su- 
perfluid and our effective theory breaks down, the classi- 
cal Boltzmann equation can be used to calculate the oscil- 
lation frequency and damping rate, as shown in Ref. |25| . 
As supplement to our main result above, we use the equa- 
tion in Ref. [2^ to calculate the oscillation frequency. As 
mentioned there, the formulas did not take into account 
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FIG. 4. (Color online) Damping rate of collective excitations 
by interpolating results from our effective field theory in low 
temperature and the Boltzmann equation in high tempera- 
ture. Red crosses are the prediction from our calculation in 
low temperature superfluid regime, and green crosses are from 
the Boltzmann equation. The blue solid line is the interpola- 
tion. The black squares are the experimental data of damp- 
ing [H. (a) l/kpa = -0.45; (b) l/kpa = -0.55. The first 
peak moves toward higher temperature when the system gets 
closer to the resonance (i.e., smaller |f/fcFt[|). 



the effect of Pauli blocking, which means that the clas- 
sical Boltzmann results are more valid at high temper- 
atures. The calculation takes into account the trapping 
potential. At 1/kpa = —0.45, if we interpolate our the- 
ory in the low temperature superfluid regime with the re- 
sults from Boltzmann equation in the high temperature 
Fermi liquid regime, two peaks by two different methods 
appear in Fig. 21 which agrees to the experimental results 
fairly well (cf. Fig. 2(b) in Ref. [ll|). Thus, we conclude 
that the first sharp peak observed in Ref. 



ll| is due to 



the superfluid phonon and fermionic quasi-particle inter- 
action, mostly through channel a in Fig. ^a) . The sec- 
ond broad peak is given by the Boltzmann equation from 
Ref. [25[ , which signals the hydrodynamic to collisionless 
transition. Our model further shows that the first damp- 
ing peak moves toward higher temperature when the sys- 
tem gets closer to resonance (i.e., smaller \l/kFa\), as 
one can see in the change from Fig. lU^b) to Fig. SJ^a). 
This phenomenon was first reported in the experiments 
of Ref. [11| (cf. Fig. 3 therein). In summary, our theory 
provides a consistent explanation for the experiments of 
damping. 

According to our calculation, a damping due to 
phonon- fermion interaction should happen in the unitary 
limit. While this was not reported in the experiments of 
Ref. [11[ , in the experiments of Ref. Q , the authors found 
that the damping rate of a Fermi gas at unitarity displays 
a weak peak immediately followed by a notch near tran- 
sition as temperature increases. Such a damping notch 



is consistent with the dip of Fig.|4]that we propose here. 

There are several reasons for the quantitative discrep- 
ancy between our calculation and the experiments. The 
most important thing is that our calculation is based on 
mean field results of Aq and yU from solving Eq. (|Iip . 
Quantum fluctuations intend to destroy the superfluid 
phase, i.e., reduce the transition temperature Tc below 
the mean field results. More reliable inputs for Ap and fi 
from calculation including quantum fluctuations [l9| will 
give a lower T^ which will make our results better agree 
with experiments. Secondly, our calculation is based on 
Fermi gas in free space while in the experiments, the 
gas is always trapped. Thirdly, we used approximation 
to solve the damping. However, to be strict one needs to 
exactly solve for the poles. Also, wq is read from Ref. [11| 
and treated as a constant in our work. That is just an ap- 
proximation since uq is also changing slightly with tem- 
perature. 

To verify our main conclusion that channel a is the 
dominating process, we propose an experiment, which is 
to measure the damping rate 7 while varying the oscilla- 
tion frequency wq, and observe how 7/a;Q changes. This 
can be done by either increasing the particle number or 
reducing the trapping potential. If channel 6, the pair 
breaking mechanism, is dominating, one should observe 
that the damping peak becomes increasingly narrower, 
and eventually is too narrow to be observable, as wq de- 
creases. The reason is that the parameter regime to sat- 
isfy the energy conservation loq ~ 2A diminishes with de- 
creasing wq, and eventually vanishes. On the other hand, 
if channel a, the Landau damping mechanism, is domi- 
nating, no matter how small ujq is, as long as Eq. (|36p is 
satisfied, the energy conservation law is always satisfied. 
Eventually 7/a;o will be independent of uiq in the limit of 
vanishing wq as discussed before in Eq. (|¥T|) . and there- 
fore the peak of damping 7/^0 remains unchanged and 
should always be observable. 

To summarize, we have studied the damping of col- 
lective modes for ultracold Fermi gases by effective field 
theory, and compared our results with experiments. By 
varying the temperature [11| our theory predicts that one 
should find two peaks of damping. The first is a sharp 
peak due to the finite-temperature phonon-fermionic 
quasi-particle interaction, and the system changes from 
hydrodynamic superfiuid to hydrodynamic Fermi liquid. 
The second is a broad peak due to the transition from 
the (coUisional) hydrodynamic to collisionless Fermi liq- 
uid (sec Fig. 4 in Ref. [l3). In the experiment that 
varies the magnetic field [8|, Il2l | but keeps temperature 
sufficiently low, one should see only the peak due to the 
phonon-fermionic quasi-particle interaction. After pass- 
ing this peak, the system is already in the collisionless 
regime without the need of going across a coUisionally 
hydrodynamic regime. Thus, there is no second broad 
peak. 
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